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ABSTRACT 



Deep dielectric charging is the suspected mechanism for 
formation of potential barriers aboard the ISEE 1 
spacecraft. Energetic electron distribution functions in 
the plasmasheet were examined for both surface and deep 
dielectric charging. Surface charging was found to be 
dependent on whether the satellite surfaces were in shadow. 
The surface potential is regulated by photoelectric 
emission, and is two orders of magnitude higher than other 
mechanisms . Deep dielectric charging deposits charge within 
dielectrics, and is independent of surface effects, such as 
photoemission and radiation-induced conductivity. 

Deposition of electrons into solar array cover cells begins 
at approximately 10 keV. 
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INTRODUCTION 



I . 



This thesis will examine spacecraft charging and, in 
particular, deep dielectric charging. Spacecraft charging 
is an important consideration in operating satellites in a 
plasma environment. A build-up of charge on the surface and 
within dielectric materials is identified with abnormal 
system behavior. These anomalies have been linked with 
unexpected system resets, false signals, and even total 
satellite failure . 

Theoretical calculations of energetic electron fluxes 
will be compared to observed in-situ data, to determine if 
charged particles present have sufficient energy to cause 
deep dielectric charging. The particle energy level for 
which differential charging begins will also be 
investigated . 

In subsequent chapters, the structure and 
characteristics of the magnetosphere will be explored. 
Several different electron distribution functions are used 
to describe the ambient plasma and how it effects surface 
and deep dielectric charging of satellites. Radiation- 
induced conductivity can alter the characteristics of 
exposed dielectrics, which in turn, can effect the re- 
distribution of deposited charge. Contributions to the 
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charging process by different mechanisms vary by several 
orders of magnitude. Comparisons will show which are 
important, and those that can be ignored. Finally, a 
correlation will be made between the calculated integral 
flux and data from geosynchronous satellite charging events. 
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II. BACKGROUND 



A. Environment 

The earth is subject to a magnetized flow of plasma 
called the solar wind (Figure 1) . Interruption of the solar 
wind flow around the earth, by the earth' s magnetic field 
forms a structure called the magnetosphere. Within the 
magnetosphere are several components which comprise its 




[Ref 1] 

infrastructure. The magnetopause separates the earth's 
magnetic field from the solar wind. Inside the 
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magnetosphere is the plasmasheet. Inward of the plasmasheet 
is the plasmasphere. The plasmasphere exists as a non- 
uniform dipolar region, centered on the earth's north-south 
axis . 

Satellites pass through three main regions daily while 
in geosynchronous orbit. On the day side, from 1200 to 1600 
local time (LT) , the satellite is generally between the 
plasmasphere and the magnetosheath, in a region populated by 
the remnants of convected plasmasheet electrons, and ion and 
electron flows out of the ionosphere (re-filling) . In this 
region, the spacecraft will transit a cool, dense plasma, 
and be continually in sunlight. At 1600 (LT) the satellite 
will move into the plasmapause, the region separating the 
plasmasphere from the plasmasheet, or on very quiet days, 
into the dusk bulge of the plasmasphere. The plasmasheet 
can inject energetic particles into this boundary region, 
raising the overall temperature of the plasma. By 2000 LT 
the spacecraft typically enters the plasmasheet, an area of 
energetic particles, and remains in this region until 0800 
LT. The hot plasma encountered (several KeV) encourages 
vehicle charging. One of the results of the charging is 
anomalies. Anomalies are erroneous satellite functions that 
are generated by the re-distribution of charge deposited by 
the plasma. These system upsets cause a wide range of 
electrical problems, including, in severe cases, satellite 
failure. Anomalies will be discussed in greater detail 
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later in this section. By 0800 LT, the electron temperature 
has dropped, as the plasmasheet electrons drift sunward. 

The probability of satellite charging drops rapidly as the 
satellite moves into the day side. Table 1 summarizes the 
composition and characteristics of the magnetospheric 



TABLE 1 



Region 


t’mJLn 


^max 


e / cm 3 


Energies (eV) 


Plasmasphere 


.2 


4 


10-10 3 


.1 - 1.0 


Flasmasheet 


8 


30 


. 1-10 


10 3 - 4x1 0 4 



plasma [Ref 1, p. 2-14 to 2-35, Ref 2, p. 3681 to 3691, Ref 
3, p. 7 to 11, p. 81 to 85] . 

The magnetosphere is populated with charged particles 
via diffusion/acceleration from the ionosphere and perhaps 
the solar wind. At quiet times, the outer plasmasphere can 
move out to include geosynchronous orbit [Ref 4, p. 3587 to 
3611] . During magnetic substorms, injections of 10 KeV to 
several MeV electrons from the plasmasheet are observed at 
geosynchronous altitude. It is thought the infusion of 
these energetic electrons is responsible for spacecraft 
charging [Ref 5, p. 278 to 308, Ref 6, p. 4-1 to 4-16]. 
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Figure 2 summarizes the movement of the regions at 
geosynchronous altitude . 




Figure 2 

B Spacecraft Charging 

Spacecraft charging is a phenomena observed when a 
spacecraft interacts with any plasma. The result of this 
interaction is an accumulation of charge within the body of 
the satellite as well as a buildup on the surface. The 
accumulation of charge on a spacecraft appears in two 
distinct ways. Absolute charging is where the potential of 
the entire spacecraft changes with respect to the ambient 
plasma. The uniform change in potential, with respect to 
the plasma will not, in general, effect the onboard 
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electronic systems. Differential charging is where 
dielectric surfaces of the spacecraft are charged to 
different potentials relative to each other. A significant 
result of differential charging is strong local electric 
fields on the surface and within the body of the satellite. 
These electric fields give rise to the possibility of 
damaging electrical discharges [Ref 1, p. 2-1 to 2-3] . 

1 . Surface Charging 

Surface charging is the buildup of charge on the 
surface of a satellite from its interaction with the plasma 
environment. The phenomena of large negative potentials was 
first recognized in the early 1970's from data taken by the 
Applied Technology Satellite (ATS) 5 [Ref 7, p. 561 to 569] . 
Surface charging at geosynchronous orbit typically reaches 
several hundred volts negative in sunlight, and to more than 
-10 4 volts in eclipse. The potentials vary with the 
environment, on very short time scales (T < 1 second) in 
eclipse, longer in sunlight [Ref 8, p. 5657 to 5667, Ref 9, 
p. 263 to 276] . 

There have been historically many instances of 
satellite charging. Deforest showed a correspondence 
between injections of energetic plasma from magnetic 
substorms and large negative potentials on Applied 
Technology Satellite (ATS) 5 when the spacecraft was in 
eclipse [Ref 7, p. 651 to 659, Ref 4, p. 3587 to 3611] . ATS 
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5 was found to charge to potentials of -5 to -10 kV in 
eclipse, and later, ATS 6 was found to charge to similar 
potentials in eclipse. In Figure 3, the charging of ATS 6 
in eclipse is illustrated. Electron and ion counts at 18.4 
keV are shown, along with the potential. The injection of 
hot plasma at approximately 21:23 UT is signalled by the 
increase in the 18.4 keV electron flux. The response of the 
satellite is to charge to approximately -6 kV in less than 
one minute. This illustrates the short time scale 
associated with charging of the satellite frame (The 
satellite capacitance to space is relatively small) . 



"HASGtNr, e Cl I fSt INJECTION 

octobcb ? ws 



A. A 



v\. A ... 



X 7 









. -*/ 



!i "n 



. lUl „ .... 
’ 1 ! 1 1** U 



1 



V,. M 



li'h" 



<*l Ut 



Figure 3 
[Ref 22] 

On 8 October, 1975, ATS 6 charged to a record negative 
potential, to 19 kV in eclipse [Ref 10, p. 362 to 366]. 
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ATS 5 and 6 both charged to negative potentials in 
sunlight. This was uncommon for ATS 5, but the satellite 
did reach a potential of -400 V on one occasion [Ref 8, p. 
5657 to 5667]. ATS 6 reached -1.5 kV [Ref 10, p. 362 to 
366] . Figure 4 shows that the time scale for sunlight 
charging is much longer than for charging in eclipse. Here 
the injection of hot plasma occurs at approximately 0800, 
but the satellite potential does not peak for almost an 
hour. The subsequent decline in potential is associated 
with the cooling of the ambient plasma. 



ATS-6 POTENTIAL: CHARGING IN SUNLIGHT 



spacecraft 

POTENTIAL, 

V 




Figure 4 

[Ref 22] 

Reasoner et al. , did a statistical study of 40 days 
of ATS 6 charging in daylight. The study revealed that the 
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probability of charging was over 50% during the local 
midnight to dawn sector as shown by Figure 5 [Ref 11, p. 89 
to 101] . 



Additionally, in a survey of high energy charging events on 
the near-geosynchronous P78-2 satellite (SCATHA) , Mullen et 
al. , found a statistical distribution for charging events 
that occurred between 1900 LT and 0900 LT . The charging 
occurred from 5.5 to 7.7 P, E . It was also shown that 
spacecraft potentials correlate to electron fluxes 
with energies of 30 keV or more [Ref 12, p. 1474 to 1490]. 

The International Sun Earth Explorer (ISEE) 1 
spacecraft was built to electrostatic cleanliness 
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Figure 5 
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specifications which it was thought would prevent daylight 
charging. However, ISEE 1 did charge to significant 
negative potentials (~ -100 V) in sunlight while near 
geosynchronous orbit in March, 1978 [Ref 13, p. 5568 to 
5578] . This charging event motivates the work done here. 
Figure 6 shows a charging event on 17 March, 1979 on board 
ISEE 1 [Ref 13, p. 5568 to 5578] . Note that the time scale 
is over 4 hours. This is 4 to 5 times as long as the 
charging period that Figure 4 shows . 




DAT 76 



Figure 6 
[Ref 13] 

Away from geosynchronous orbit, satellite charging 
occurred on the Defense Meteorological Satellite Frogram 
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(DM3F) satellite while flying through auroral regions. 
Negative potentials from 50 V to 1.5 kv were observed while 
the satellite crossed high magnetic latitudes in low earth 
orbit (LEO) [Ref 14, manuscript]. 

2 . Differential Charging 

Differential Charging is the charging of a piece of 
a spacecraft relative to an adjacent section. Differential 
charging occurs, in general, because of variations in 
materials, shadowing of surfaces, and the deposition of 
charge into dielectrics . There are a variety of materials 
used for satellite surfaces, many of which are non- 
conductors (dielectrics), such as kapton, teflon, and glass. 

Photons from the sun incident on a spacecraft cause 
electrons to be emitted from the surface, provided the 
energy of each photon is greater than the work function of 
the material- Because of this shadowing is important. The 
nature of such shadowing depends on satellite geometry and 
stabilization technique- The most popular techniques are 
three-axis and spin stabilization. In the latter case, the 
orientation of the spin axis relative to the sun-satellite 
line is important. 
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As previously mentioned, surface charging in eclipse 
can reach kilovolt potentials, while the satellite potential 
in sunlight is generally not as large [Ref 8, p. 5657 to 
5667]. Shadowed dielectrics will behave as though they are 
in eclipse, and charge to much larger negative potentials 
than the satellite frame or illuminated dielectrics. A 
potential gradient of kV/cm across dielectrics can lead to 
arcs and electrical discharges if the accompanying electric 
fields exceed the dielectric breakdown strength [Ref 5, p. 
277 to 308] . 

The effects of differential charging have been 
observed- ATS 6 data revealed evidence of differential 
charging. Figure 7 shows an electron distribution observed 
by ATS 6. [Ref 15, p. 715 to 719]. This figure shows that 
electrons with energy below 50 eV must have been generated 
on or near the satellite surface. Hence, the photoelectrons 
and secondary electrons generated on the satellite surface 
were being reflected by some form of barrier, and were 
unable to escape the vehicle. Otherwise, astronomical 
values for the plasma density result. Olsen et al . , 
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postulated the barrier was caused by differential charging 
vice electron space charge [Ref 16, p. 6809 to 6819]. 




Figure 7 



Differential charging of dielectrics could be 
directly measured with the satellite surface potential 
monitor (SSFM) on the SCATHA satellite. Dielectric 
materials (kapton, teflon, Si0 2 ) were found to charge to 
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negative kilovolt potentials when shadowed as shown in 
Figure 8 . 
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Figure 8 
[Ref 18] 

This differential charging occurred with a time 
constant of tens of seconds, as appropriate to the 
capacitance associated with the thin dielectrics. This is 
in contrast to the more rapid charging of the satellite as a 
whole (t < 1 second) , which has a relatively lower 
capacitance to space [Ref 18, p. 506 to 509]. The 
differential charging is also responsible for the charging 
of the vehicle as a whole [Ref 17, p. 14-1 to 14-12]. 
Modelling of daylight charging showed that the barrier 
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induced by differential charging allowed a net negative 
current to the sunlit surfaces. 

C. Anomalies 

The magnetosphere is often an unfriendly environment in 
which satellites must survive and operate. Spacecraft are 
subjected to high energy particle fluxes, and discharges 
resulting from surface and differential charging can produce 
electromagnetic pulses (EMP) . These EMP can enter into the 
internal circuitry through unshielded cables, ungrounded 
components, and any external conduit that will transport the 
Fourier spectra of the EMP. Additionally, very high energy 
particles can cause single event upsets (SEU) [Ref 1, p. 4- 
22 to 4-23] . 

Satellite anomalies were attributed to charging 
following the realization in the mid-1970's that large 
negative potentials could be found on satellites positioned 
at 6.6 earth radii. Figure 9 shows a distribution of 
anomalies as a function of time of day for several 
geostationary spacecraft. This historical figure was 
promptly related to the charging distribution obtained from 
ATS 6 (see Figure 5) and a correspondence was suggested [Ref 
11, p. 89 to 101] . 

A variety of anomalies for several different satellites 
are tabulated in Reference 1. They include power-on resets, 
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EMP, spin-up, false flags in software, clock shifts, false 
or phantom commands, telemetry errors, power-downs, 
circuit noise, power system shutdowns, system upsets, and 
total satellite failures. The figure is a distribution of 
anomalies observed on several geosynchronous satellites. 
The times are local and the distance from the origin is 
solely to allow visual resolution of the events [Ref 1, p. 
5-2 to 5-11] . 



Figure 9 
[Ref 11] 

Although there are many examples of anomalies in 
geosynchronous satellites, most are difficult to relate 
unambiguously to satellite charging as a cause. The most 
revealing case was with the SCATHA satellite on 22 
September, 1982, when differential charging, EMPs, and 
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anomalies were found together as shown in Figure 10 [Ref 19, 
p. 239 to 243] . 




a) The absolute value of the voltage on the Satellite Surface 
Potential Monitor gold sample as a function of time. The potential of 
the sample Is negative with respect to the ground reference; b) the 
count rate for 18.4 keV electrons as a function of time. 




UNIVERSAL TIME, hrimin^ec 



The timing sequence of the antenna switch on the VLF Analyzer 
experiment during one time period when anomalous behavior o ccu rr e d on 
September 22, 1982. 



Figure 10 

During a period of energetic electron flux and high 
surface potentials, the timing sequence of a switch in an 
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experimental package showed anomalous behavior. The arrows 
at the bottom of the figure show a flip in logic states 
coinciding with discharge events. 

An example from a completely different environment, the 
Jovian magnetosphere, provides an unusually clear 
demonstration of the relationship between anomalies and an 
environmental cause. A sequence of anomalies occurred when 
Voyager I passed Jupiter on 5 March, 1979, as seen in Figure 
11 . 
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Figure 11 

[Ref 20] 

Prior to and during the fly-by at 5 R Jf a series of 
power-on resets (POR) plagued a critical photographic 
sequence. Figure 11 shows resets and the energetic electron 
flux. There is a clear correspondence [Ref 20, p. 21 to 
30] . 
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Laboratory simulation of the POR process identified 
unterminated spare leads in a cable bundle in which the POR 
cable co-existed. A portion of this cable exited the 
spacecraft body and was exposed to the plasma environment. 
Re-creation of the observed environment produced 25 mA 
current pulses capable of inducing the POR. Surface 
charging in a cold dense plasma was ruled out and a 
correspondence was established between the POR and the 10 
MeV e~ flux component measured by on-board instruments [Ref 
20 , p. 21 to 30] . 

This example shows the importance of energetic 
electrons, the universality of spacecraft charging, and 
anomalies. It is also one of the prime examples of deep 
dielectric charging, which is the deposition of charge into 
dielectrics (like cable insulation) where the charge is 
trapped, accumulates, and ultimately is released as an EMP . 
Deep dielectric charging is the primary focus of this 
thesis . 

Reagan et al . , were among the first to discuss deep 
dielectric charging (Figure 12) . They presented SCATHA data 
for days with large fluxes of high energy electrons (several 
MeV) and found a correspondence to charging. Energetic 
electrons incident on satellite surfaces can bury themselves 
either into dielectric materials or into floating metal 
(conductors not connected to the spacecraft ground) . As 
charge builds up within the satellite body, it may exceed 
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the breakdown strength of the material and arc. The 
associated EMP may be sufficiently large to induce anomalous 
behavior in one or more of the satellite electronic 
components . 



ELECTRON FLUX 




CHARGE TRAPPED 
ON ISOLATED 
CONDUCTOR 




Figure 12 
[Ref 1] 

The deposition and buildup of charge in the dielectfic 
were modelled and found in agreement with in situ 
measurements. Although peak electric field strength was 
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below material breakdown levels, rapid changes in these 
fields resulted in large EMP [Ref 21, p. 354 to 365] . 

Baker et al . , have also looked at the consequences of 
high energy electron fluxes. Figure 13 shows a connection 
between high energy electrons at geosynchronous orbit 
(measured on 1982-019) and failure of a Geosynchronous Orbit 
Earth Satellite (GOES) 5 encoder lamp. 




Figure 13 
[Ref 6] 

Another clear indication of the high energy electron 
flux-anomaly cause and effect relationship is shown in 
Figure 14. Over a period of two years, spacecraft 1979-053 
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experienced multiple malfunctions in its star tracker 
system, as indicated by arrows above 3 MeV fluxes. Notably , 
anomalies occurred only when the electron flux exceeded 
approximately 6 count rate units. After several days of 
erratic upsets, the anomalous behavior disappeared [Ref 6, 
p. 4-1 to 4-16] . 




Figure 14 
[Ref 6] 



D . Summary 

Surface and differential charging are common occurrences 
on geosynchronous satellites. Deep dielectric charging can 
also occur in the presence of high energy electron fluxes. 
These mechanisms are important aspects of spacecraft 
operation and survivability. This thesis will build on the 
work done in previous examinations of spacecraft charging 
and extend it . 
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Ill . 



THEORY OF CHARGING 



A . The Environment 

Charging is the buildup of charge on the surface and 
within the body of a spacecraft immersed in a plasma. The 
plasma environment considered here is the magnetosphere , 
particularly at geosynchronous orbit. The plasma at 
geosynchronous orbit has it's origins in the solar wind and 
the earth's ionosphere. These sources of electrons and ions 
continually introduce energetic particles into the earth' s 
geosynchronous altitude region where stationary satellites 
orbit - The energies of the particles vary from fractions of 
an ev to MeV. Magnetospheric storms and sub-storms inject 
high energy particles into this region. This is especially 
common from 2000 LT to 0900 LT, when geosynchronous 
spacecraft transit the plasmasheet . 

The distribution function, or phase space density, 
characterizes a plasma. Equation (1) shows a Maxwellian 
distribution, as a function of the velocity of the 
particles. Such distributions are often used to describe 
the plasma observed at geosynchronous orbit. 
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( 1 ) 



f ( v) = n 



1 77 

2nkT 



, 3/2 




The temperature, kT, is a measure of the average energy of 
the distribution, where k is Boltzmann's constant, and T is 
the temperature in Kelvin. Electron volts (eV) are also 
commonly used units. The mass of the species is m, the 
velocity of the particles is v, and the density n. 

Flux, or flow, denoted by J 0 , is equal to the product of 
the particle charge, density, and velocity. We will need to 
relate J to the currents and the particle detector 
measurements- The concept of flux can be further defined by 
including the originating direction of the particles. Then 
the equation j 0 = J 0 /4 K defines a portion of the omni- 
directional flux arriving from a particular direction. This 
implies that the surface is always perpendicular to the 
incident flux. In order to account for the flux from 
directions that are not normal to the surface, the equation 
must be integrated over all angles from which particles are 
incident onto the surface, denoted by dA, where dA is an 
outward vector normal to the surface. The particle 
detectors used on spacecraft, typically electrostatic 
analyzers, are designed to measure a wide range of energies 
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by sequentially sampling a specific energy range over a 
small angle of incidence. Consider one energy channel to be 
between and energy E and E + dE, along a solid angle, 
denoted by d£l, and over a time period from t to t + dt . 

These parameters will then define the flux of particles 
within a specified energy range, from a specific direction, 
at one location, over a distinct time interval incident onto 
a surface area of the satellite. The measurement in each 
energy channel gives differential energy flux j [Ref 3, p. 

36 to 37]. Figure 15 illustrates the angles associated with 
differential flux measurements. 



Z 




Figure 15 
[Ref 3] 

Equation (2) expresses the differential flux. 
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( 2 ) 



dN - j ■ dA dE dQdt 



The units of differential flux are # / (cm 2 s ster KeV) . 

To find the distribution function from the differential 
flux, consider non-relativistic particles in velocity space. 
The relationship between energy and velocity is given by the 
equation v 2 = 2-E/m. Using this relationship, the number of 
particles in a solid angle dQ in the energy range dE can be 
expressed in a velocity space volume by v 2 dv dii. The 
number of particles in a unit volume of space is then given 
by the distribution function f multiplied by the number of 
particles with velocity v in a time t incident on an area 
dA. Since this is the same as equation (2), the 
differential flux j is therefore related to the distribution 
function f by the equation j = f -v 2 /m [Ref 3, p. 36 to 38] . 

This thesis assumes that space is isotropic, that the 
distribution of particles does not have a preferred 
direction. By assuming an omni-directional flux, the 
distribution of particles can be characterized by a single 
Maxwellian, where particle velocities are independent of 
orientation. This characterization of the particle 
distribution gives an accurate description of the plasma at 
lower energies. When the particle flux has a high energy 
component out on the tail of the distribution, other 
mathematical descriptions are used in conjunction with the 
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Maxwellian. The higher energy particles used herein are 
characterized by an inverse power law. 

Further assumptions about the plasma are that it exists 
in a charge-balanced state. This concept of quasi- 
neutrality states that the overall charge of the plasma will 
be neutral. Additionally, the Debye length is large enough 
for statistical concepts to have meaning. This implies that 
collisions are rare and particle interaction can be 
neglected [Ref 3, p. 22] . 

From kinetic theory, the average energy of a species in 
an ensemble can be represented by 



( 3 ) 



Ei = — in.-vj 

i 2 1 




where Ei is the average kinetic energy, v A is the average 
velocity, and T ± is the temperature, all of the ith species 
[Ref 3, p. 27]. This implies that 



( 4 ) 
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KT 

m 



From equations (3) and (4) , it can be seen that in an 
equi-temperature species, the electron velocity is at least 
43 times the ion velocity, by virtue of their light mass. 
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Therefore, for comparable temperature, the effects of ions 
can largely be ignored. 

B . Surface Charging 

Spacecraft are constructed from conducting and 
insulating materials. Metal surfaces provide a ready 
pathway for incident charges to re-distribute themselves. 
Insulators or dielectrics, on the other hand, do not easily 
allow charges to move within or over them. As particles of 
varying energy hit the satellite surfaces, some strike 
conductors and move to establish equi-potential surfaces. 
Other particles collide with dielectrics and are fixed, thus 
establishing electric fields over localized portions of the 
spacecraft surface. These electric fields are dynamic in 
nature, always trying to move the charge to positions of 
minimum energy. 

The distribution of charge and currents on the 
spacecraft surface at equilibrium will minimize the 
potential energy. This implies that any charge on the 
satellite surface will, if mobile, arrange itself, with 
respect to the plasma currents and surface potential, at 
minimum energy . 
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The balance of currents entering and leaving the surface 
of a satellite can be characterized by equation (5) 

L ^ambient + ^ energetic^ ^ secondary* ^backseat ter + m 

True ambient flux is all of the flux that surrounds the 
spacecraft. Ambient current density is current/unit area 
that comes from the ambient flux. 

A distinction is made here regarding the ambient and 
energetic fluxes. Even though both come from the overall 
ambient flux, there is an artificial break point established 
at 40 keV. At this point, the description of the 
distribution function changes from a Maxwellian to an 
inverse power law. Therefore, it is easy to separate the 
two current densities in a similar way. 

The secondary current density comes from incident 
electrons that knock lose one or more electrons from the 
material surface. Electrons that are reflected from the 
surface are included in the backscatter term. Finally, the 
photoelectric term accounts for electron emission of bound 
electrons from the surface (photon energy > material work 
function) . 

Figure 16 shows how each of the current densities in 



30 



equation (5) contributes to the net current density at the 
spacecraft surface. 



PHOTON 




Figure 16 
[Ref 23] 

Equation (6) shows that the current density is 
proportional to both the electron density and the 
temperature . 



< 6 > J = g n v chermal 

where q is the charge of the species, n is the density, in 
particles / volume, and V thsrTnal is given in equation (4) . 
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The net current density is dependent on many different 
parameters. The temperature of the ambient and energetic 
distributions,, plasma density, velocity of the spacecraft in 
the plasma, percentage of projected area that is dielectric, 
materials from which the dielectric surfaces are made, and 
whether the satellite is in sunlight or eclipse, are all 
factors in determining J net . 

Additional sources of current come from secondary ion 
emission, movement of the satellite across local magnetic 
fields or relative to the ambient ion population, ions from 
satellite outgassing, and induced current flows from high 
differential potentials. These are usually minor effects 
that will be ignored. [Ref 1, p. 2-1, Ref 24, p. 168] 

Equation (5) describes a delicate balance between 
different mechanisms that contribute to the surface charge 
of a spacecraft. In turn, each of the contributions to J net 
will be examined. 

1 . Ambient Current Density 

Current density equals the charge multiplied by the 
flux. The flux is found when differential energy flux is 
integrated over the applicable energy range. The flux is 
given by 

= J'JJ'f(v) ■ ( v • n ) d 3 v 



( 7 ) 



where the f (v) is given by 



(i) 
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The particles incident on the surface area from non- 
normal angles is accounted for in the (v • n) term in 
equation (7) [Ref 22]. Using spherical coordinates, the 
flux in velocity space can be expressed as 



( 8 ) 
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V 3 COS 6 sin <9 d(f) d 6 dv. 



The angular integration gives a factor of 71. By using the 
relationship E = ^mv 2 r the integration can be done over the 
distribution function energy range. This gives 



( 9 ) 
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The flux is the integral of the differential energy flux, 
J (E) , integrated over the energy range. The differential 
energy flux is 
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(10) 



J{E) = n 



2ti m 



— \ 1/2 e 




2. Secondary Emission Current. Density 

Electrons incident on a spacecraft surface can 
undergo anyone of three interactions, reflection, 
backscattering, and secondary emission. Reflection of 
electrons is a extremely low energy phenomena, generally 
found at energies below 10 eV, and can be neglected [Ref 
23] . 



When low energy electrons collide with the satellite 



surface, they quickly come to rest. The energy from the 
incident electrons becomes available and may be absorbed by 
bound electrons. If this energy is sufficiently high, one 
or more electrons may be ejected from the surface material. 



Secondary electrons will be ejected at a low energy, 



generally with an energy below 50 eV. More than 80% of all 
secondary electrons have energies below 20 ev [Ref 30, p. 
163 to 189] . A secondary yield term, attributed to 
Sternglass (1950), is material and ambient flux dependent. 



( 11 ) Yield s<iconds[v = 7.4 -6 



max 
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is the greatest 



where E ma „ is the energy where yield S9Condary 
and is the secondary yield at E max . Table 2 shows 

secondary emission properties for various spacecraft 
materials [Ref 23, p. 41, Ref 25, p. 123, Ref 26, 493 to 
499, Ref 27] . 



TABLE 2 



MATERIAL 


E (max) eV 


6 (max) 


Gold 


800 


1 . 45 


Oxidized Aluminum 


300 


2.60 


Silicon Dioxide 


420 


2.50 


Aquadag 


350 


0 . 75 


Teflon 


200 


3.00 


Kapton 


250 


1.80 



Figure 17 shows a graph of the secondary electron 
yield for kapton as a function of the incident electron 
energy . 
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The figure is typical for most spacecraft materials. 
The maximum yield on the curve, usually from an incident 
electron at a few hundred eV, results in more than one 
electron ejected from the surface. 

In the figure, secondary emission electrons 
outnumber incident electrons between 50 eV and 1 keV. These 
two energies, where the secondary yield is one, are known as 
the lower and upper crossover energies, respectively. If 
the ambient flux is sufficiently rich in particles between 
the two crossover energies, more electrons will leave the 
surface than are deposited, J secondary > Jambi.nt' and the 
spacecraft will charge positively, typically a few volts 
positive . 

The upper crossover energy suggests the threshold 
energy effect- When the electron distribution function 
contains sufficient high energy components and of sufficient 
energy to balance the out-flux of secondary electrons, the 
satellite potential will become negative. Therefore, the 
temperature of the ambient electron distribution will be 
instrumental in the determination of the spacecraft 
potential. The energy at which balancing occurs is known as 
the threshold energy and is the start of negative charging. 

The threshold energy will vary among spacecraft 
because of geometry and dielectric materials, but lies 
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SECCNCA3Y EMISSION YIELD (5) 



within the 10 to 20 keV band [Ref 26, p. 493 to 499, Ref 39, 
p. 169 to 175] . 




rniMAnY EitcinoN eneiioy |pvi 



Fi.gure 17 
[Ref 26] 
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3 . Backscattered Current Density 



The ambient flux electrons have different kinetic 
energies. As electrons hit the satellite surface, some 
reflect, others cause secondary emission, still others are 
backscattered. Backscattering is primarily elastic 
collisions (coulomb scattering) of the incident particles 
with the surface atomic nuclei. Low-energy electrons are 
deflected into a wide scattering cross-section. The 
responsible mechanism is complicated, but a useful model for 
backscattered electrons is [Ref 25, p. 104 to 105, Ref 29, 
p. 1969 to 1981] 

(12) Yiel d backscattered = 0.1 (l.O - ( P + e ) 



P is given by 



( 13 ) 



P = 10.0 



1.0 - ( AiO ) <°- 037 *> 



where E is energy in keV and Z is the effective atomic 
number. This formula is good for incident electrons with 
energy > 50 eV. Incident electrons with less energy are 
counted along with secondary emissions. Figure 18 shows how 
the backscatter emission curve determines the backscatter 
flux for a Maxwellian distribution with T = 7.7 keV and n = 
0.15 cm' 3 . 
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Energy (eV) 



Figure 18 
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Two transcription errors were noted in References 23 
and 25. The exponential term in P was not raised to a 
power, and the 0.05 was written as 0.5. This could have 
resulted in a very slight under-calculation of the 
backscattered flux between 100 and 500 eV in the referenced 
work . 

4 . Photoelectric Current Density 

When photons are incident on a spacecraft, electrons 
are ejected from the surface , provided that the energy of 
the photons is greater than the electron binding energy of 
the surface molecules . Typically, the photon energy 
spectrum is higher than the threshold energy or work 
function of the material. The photoemission electrons, 
therefore, have an energy spectrum (1 to 15 eV) that 
modulates with the photon spectrum, but is lower by the 
threshold energy. Saturation current densities range up to 
4 2 (.lA/nf 2 . 

The wavelength of the incident photons is inversely 
proportional to their energy. Figure 19 shows a comparison 
of the solar flux to photon wavelength, and the resultant 
photoemission yield. A X of 10 4 A corresponds to a photon 
energy of 1.2 eV, 3*10 3 A to 4.1 eV, and 10 3 A to 12.4 eV. 

For tungsten, incident photons > 10 eV eject an electron 1% 

~ 25% of the time [Ref 31] . 



40 



Table 3 lists saturation current densities for 
spacecraft materials , mean kinetic energy upon ejection from 



the surface, and the electron density next to the ejecting 
surface . 




Figure 19 
[Ref 30] 
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PHOTOEMISSION YIELD 



Table 3 summarizes computed photoemission values for typical 
spacecraft materials [Ref 30, p. 163 to 189]. 



TABLE 3 



Material 


Saturation J 




e" density iurf 




(|XA/m 2 ) 


(eV) 


(cm" 3 ) 


A1 2 0 3 


42 


1.33 


1670 


Indium Oxide 


30 


1.35 


1260 


Gold 


29 


1 . 40 


1110 


Aquadag 


18 


1 . 17 


500 


Cr. Steel 


20 


1 . 42 


770 


Graphite 


4 


1 .40 


150 


Average 


21 


1 . 42 


830 



5 . Energetic Current Density 

The energetic current density is that portion of the 
electron spectrum which is no longer well-described as a 
Maxwell-Bolt zmann distribution. We use the term energetic 
here to signify that the differential flux component comes 
from a power law distribution, which begins at approximately 
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40 keV for data considered here. The distribution function 



is proportional to (E/E 0 )~ x , where x typically 4 to 6. In 
our case, x = 5 (this gives a flux that in proportional to 
E -3 ) . The flux is found when differential energy flux is 
integrated over the applicable energy range. 

(7 ) <X> = J J J f { v) * ( v • h ) d 3 v . 



Now f (v) is given in terms of energy by 



(14) 
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The distribution function is given in terms of energy. The 
mass is m, q is the charge, E 0 is the minimum energy for the 
distribution function, in eV [Ref 13, p. 5568 to 5578, Ref 
22]. The normalization is such that n = /// f (E) when done 
with limits E = E 0 to infinity. 

Inserting the distribution function into equation 7, 
and using the relationship E = 3 ^mv 2 , the integration can be 
done over the distribution function energy range. The upper 
limit of integration is infinity. Since the distribution 
function tail asymptotically approaches zero, the 
contribution from the upper limit is zero. Therefore, the 
energetic flux is primarily determined by the lower limit of 
integration . 
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Inserting equation (16) into (7) gives 



( 15 ) 



$ = n 



7 
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As with the ambient differential energy flux (equation 11) , 
the energetic differential energy flux is given by 



current density (positive is defined as into the satellite 
surface) are the ambient flux, which is the Maxwellian flux 
from 0 to 40 keV, and the energetic flux, which accounts for 
the energetic electrons above 40 keV. Countering the 
positive flux are the secondary, backscatter, and the 
photoelectric current densities. 



(16 ) 




To summarize, the positive contributions to the net 
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Table 4 gives an order of magnitude comparison of 
which current densities are important and which mechanisms 
can be ignored. 



TABLE 4 



CURRENT DENSITY SOURCE 


FLUX (e/m 2 s) 


J (|0A/m 2 ) 


Ambient (<b < 4 0 keV) 


3.0 x 10 12 


.48 


Energetic (<£> > 40 keV) 


1.4 x 10 11 


CM 

O 


Backscatter 


3.2 x 10 11 


. 05 


Secondary 

• 


1.6 x 10 12 


.26 


Photoemission 


1.3 x 10 14 


21 . 



The largest contributor to the net current density 
is the photoelectric term, by approximately two orders of 
magnitude. When the spacecraft is in eclipse, however, the 
ambient flux provides a net accumulation of electrons to the 
surface, and the vehicle would charge negatively. This is 
closely followed by the secondary emission current density. 
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The distribution function used here is for a hot 



plasma distribution , where the temperature is higher than 
for a quiet day. During quiet days, when the ambient and 
energetic fluxes would be lower, the sum of the secondary 
and backscatter terms would be closer to the positive 
fluxes. This would allow the net current density to self- 
balance, precluding large negative potentials from forming. 

The effect on the potential from the net flow of 
charge to the surface of the spacecraft is clearly shown in 
Figure 20. On 20 September, 1974, ATS 5 and ATS 6 were 
immersed in the same plasma, and both charged to similar 
negative potentials, prior to the active experiments which 
were the focus of this work. An experiment on spacecraft 
potential modification was conducted with the electron 
source on ATS 5. The electrons came from a filament source, 
from - 0632 to 0635 [Ref 38, p. 527 to 532] . Prior to 0632 
and after 0635, ATS 5 and ATS 6 charged to the same 
potential . 

This shows that although surface material 
characteristics are important (e.g., secondary yield), 
variations between these vastly different vehicles did not 
lead to greatly different potentials in eclipse. Hence, 
simple assumptions about yield terms, such as assigning 
"typical 1 ' values for glass or aluminum are not unreasonable. 
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Figure 20 
[Ref 38] 
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C. Deep Dielectric Charging 

Energetic electrons lie in the high-energy tail of the 
electron distribution function. As the average temperature 
of the distribution increases, the number of high energy 
particles increases, and so does the likelihood of deep 
dielectric charging . 

In 1976, a new charging mechanism was proposed in which 
high energy electrons imbedded themselves in dielectric 
materials [Ref 36, p. 237 to 246]. This mechanism, deep 
dielectric charging, becomes important when the ambient 
electron distribution function has a high-amplitude high- 
energy tail. Through secondary and photoemission, or solely 
through secondary emission in eclipse, the surface layer can 
become electron depleted. This, in turn, creates large 
electric fields within the dielectrics as shown in Figure 
21 . 




Micrometer* 



Figure 21 
[Ref 37] 
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When the top layer is ’’hole” rich, the electric fields 
that are produced are on the order of 10 6 volts per cm. As 
time increases from minutes to hours, the potential will 
draw electrons within the dielectric toward the surface. 

The consequence of "energetic electrons” in the 
distribution function is deep dielectric charging. 

Important results of deep dielectric charging are energy 
deposition profiles , radiation- induced conductivity, 
localized formation of strong electric fields, and formation 
of potential barriers. 

The build up of charge within the dielectrics is 
governed by two competing processes, the rate which charge 
is deposited, and the rate which it leaks out. The net 
charging rate can be modelled by considering the dielectric 
to be a one-dimensional slab of thickness d. If the current 
density entering is J ± and J c is that which exits the 
dielectric, then the divergence of the current density is 
approximately (J ± - J Q ) /d. The continuity equation 

(17) -|e = V ■ J D + V ■ oE 

dt D 

can then be used to model the net rate of accumulation. 
Equation (18) shows the mechanism for the net buildup of 
charge, but does not include surface conductivity effects, 
as modified by radiation. 
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Solving the continuity equation for charge density as a 
function of time , assuming no charge at t = 0, as would be 
the case when the spacecraft environment is calm and there 
is photoemission, gives 

(18) p(t) = Ji ~ J ° - [ 1 - e""' t ] 

a a 

These simple equations show that the time constant X = 
e/a. For glass, t is many hours since deposited electrons 
will remain practically forever. 

1 . Radiation-induced Conductivity 

The conductivity in equation (18) is a free 
parameter in the charge balance equation. The conductivity 
of the dielectric can be altered over time by a phenomena 
known as radiation-induced conductivity (RIC) . This effect 
varies with the depth of the material, and is confined 
primarily to the surface of the dielectric. The effects of 
RIC will be discussed below. First the electron deposition 
process must be considered. 

2. Energy Deposition Profile 

The energy deposition process is a complex mixture 
of elastic and inelastic collisions, scattering, and 
material penetration, as described by scattering theory [Ref 
32, p. 2-1, 2-2]. Figure 22 shows the deposition depth of 

electrons in silicon glass as a function of the incident 
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particle energy [Ref 35, p. 123] - The significance is that 
for a 12 mil (0.03 cm) cover cell, electrons from several 
hundred eV to approximately 225 keV will deposit themselves 
within the dielectric. Electrons with higher energies will 
pass through and enter the spacecraft body. At the low 
energy limit (10 keV) electrons reach approximately 1.5 x 
10~ 4 cm, or 1.5 [A into the material. We take this as being 
the minimum energy limit for charge deposition. Electrons 
deposited at lesser depths will flow rapidly along the 
surface . 
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Figure 22 
[Ref 35] 
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The effective conductivity of the dielectric is a 
combination of the intrinsic conductivity <J L and the 
radiation- induced conductivity <J RIC , and can be expressed as 



( 19 ) 



® effective ®i + ® RIC 



[Ref 33, p. 2282]. Radiation-induced conductivity can be 
modelled by 



( 20 ) 



o 



RIC 



K D A 



where K is the coefficient of RIC, and a is a material 
dependent parameter. Although theory predicts that both are 
constants and are independent of type and energy of the 
radiation, there is disagreement on their empirical values. 
a is frequently found to lie within 0.9 to 1.1, and can be 
taken as 1.0. K lies within a greater range of minimum and 
maximum values, hence introduces greater uncertainty into 
equation (20) than does a. 
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Table 5 lists maximum and minimum values for typical 
spacecraft dielectrics. Units are 10" 19 sec/ (fi cm rad) [Ref 
34, p. 580] . 



TABLE 5 



Material 






Kapton 


60 


1.2 

I 


Teflon 


1000 


20 


Mylar 


2 . 1 


1 . 8 



The dose rate in equation (20) is given by dD/dt, 
and is the rate at which energy is deposited into the 
material, in units of rads per second. A rad is defined as 
100 ergs per gram. The dose rate consists of three parts, 
the collision stopping power, the specific thickness, and 
incident flux. 

The collision stopping power (CSP) is a measure of 
how far into a material electrons will penetrate as a 
function of their incident energy and the thickness of the 
material. By dividing the CSP by the density of the 
material, the CSP becomes independent of the material and 
comparative analysis can be performed. At low energies, the 
CSP measure loses its usefulness, and ceases to be 
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Collision Stopping Power (MeV cm /gm) 



meaningful at energies below several hundred eV [Ref 35, p. 
20] . The units of CSP are MeV-cm 2 / gm. Figure 23 shows 
CSP as a function of the energy for borosilicate, "pyrex" 
corning 7740 glass, density = 2.23 gm/cm 3 . This is similar 
to the material used for solar array cover cells. 

100 
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0 0.05 0.1 0.15 0.2 0.25 

Energy (MeV) 

Figure 23 
[Ref 35] 

The dose profile in a low atomic number (low-Z) 
material is directly proportional to incident kinetic 
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energy. By defining the incident thickness in terms of e" 
/cm 2 -MeV, the dose profile will be independent of the 
particle's incident energy and the nature of the material. 
This thickness is defined as specific thickness (T s ) r and is 
given by 



( 21 ) 




where f A is the weight fraction of the ith element of the 
material, and is given by 



( 22 ) 




A i 

Kn i Aj: 



and where the parameters are defined as: 

- p = material density (gm/cm 3 ) 

- x = depth in the material in cm 

• N 0 = Avogadro number 

• T = e” energy in MeV 

- A = Atomic weight 

• n i = number of atoms in the ith element of the material 

The incident differential energy flux is calculated from the 
plasma parameters in units of eT / (cm 2 s ster MeV) . 
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The front surface dose is found by multiplying the 
differential energy flux with the CSP . The dose profile in 
a low Z material is calculated by finding the front surface 
dose and then modulating it with the curve in Figure 24. 

The dose rate is then calculated by integrating the dose 
profile over the applicable energy range. 

Figure 24 plots the incoming flux, or relative dose, 
against specific thickness (T.) . 




Figure 24 
[Ref 32] 
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If the saturation electric field at any time exceeds 
the dielectric breakdown threshold, as modified by RIC, 
charge re-distribution or arcing will occur. 

Using a teflon sample (d = 10 mils) , a breakdown 
threshold of 10 7 coulombs/cm 2 and a discharge time of 
roughly 50 nanoseconds, peak structural currents (geometry 
dependent) can range from 1 to 10 amps [Ref 33, p. 2281 to 
2284] . The use of teflon, a dielectric material, is for an 
order of magnitude picture of the peak currents that can 
occur . 

D . Potential 

The overall charge of the spacecraft determines the 
potential. The charge density on the spacecraft is 
determined by the net current density, as shown in equation 
(5), integrated over the appropriate time interval. We 
tacitly assume the currents are constant over the periods of 
integration . 



( 5 ) 



^ne t ^ ambi ent + ^ energe tic + ^ secondary + ^ ba cksca tt er + ^hv 



For the two different charging mechanisms identified, 
surface and deep dielectric charging, the appropriate time 
constants are very different. 
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For surface charging, the time periods which are 
appropriate are typically a second or less. This is because 
the relatively small capacitance associated with satellite 
potentials can lead to significant charge buildup (changes 
in potential) within a millisecond time frame. If the 
satellite is near equilibrium, the net current will be 
small, and it should be possible to estimate the net charge 
density accumulated in one spin period by taking the 
currents to be constant over that period. 

( 23 ) Q/m 2 sul . f = (J amb - J photo /2 - J sec - J bs ) ■ spin time 



Sunlight shines on approximately 20% to 25% of the 
spacecraft surface at any given time, depending on the 
vehicle geometry. The spin time for a satellite can range 
from less than a second to about one minute per revolution. 

In deep dielectric charging, the charge density comes 
from the energetic current density. Therefore, the 
deposition of charge is dependent on the energetic tail of 
the distribution function. This gives a net deposition of 
charge into the dielectric material over long periods of 
time, that is minutes to hours. 

The deposition of charge into the dielectrics will occur 
slowly, but can eventually become significant if the 
resistivity is high. Hence the charge distribution can be 
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similar to that induced with shadowed differentially charged 
surfaces [Ref 40, p. 6809 to 6819] . 

A potential barrier around the spacecraft can form and 
grow large enough to suppress low-energy secondary and 
photoemission electrons. These electrons are returned to 
the spacecraft and help maintain a negative potential, even 
in the presence of sunlight. This is the process generally 
invoked for daylight negative potentials for satellites with 
insulating surfaces [Ref 28]. This development of 
differential potentials, or the buildup of trapped charge is 
a slow process, usually on the order of five to ten, or 
more, minutes [Ref 8, p. 5657 to 5667]. It is possible that 
deep dielectric charging could also lead to barriers, and 
also provide a mechanism for daylight charging [Ref 23] . 

This is what we want to study further. 

If we take the potential of a spacecraft to initially be 
zero, then in one spin of the vehicle the surface charge 
density, G photo , develops from incident sunlight. Can enough 
charge be trapped so that 



If so, we can imagine that a surface starting at <}> “ 0 will 
at some point develop a net charge which can not be balanced 
by photoemission in one spin period. 



( 24 ) 




^ ® photo ? 
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E . Summary 

Lower energy ambient flux is responsible for surface 
charging. This phenomena occurs in eclipse, when the 
photoelectric emission is blocked. Surface charging appears 
in a very short time span (x < 1 second), charges uniformly 
to several keV, and returns rapidly to equilibrium potential 
when photoemission is restored. 

A second mechanism, deep dielectric charging, occurs 
when there is an energetic tail on the distribution 
function. These energetic electrons bury themselves into 
the dielectrics covering the satellite surface. Movement of 
the electrons out of the dielectrics is slow (x can range 
from seconds to several minutes) . If a buildup of charge 
exceeds the breakdown strength of the material, electrical 
discharges may cause anomalies in satellite systems. 

Associated with deep dielectric charging is the 
formation of potential barriers. Upon entering sunlight, 
photoelectrons emitted from the satellite surface are not 
energetic enough to overcome the potential barrier, and are 
returned to the satellite, thus maintaining the negative 
charge. Negative potentials in the range of several hundred 
volts will remain until the charge buried within the 
dielectric can leak out. 
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IV. RESULTS 



A. ISEE 1-Observed Integral Flux 

1 . Ambient Differential Energy Flux 

On 17 March, 1978 (day 76), ISEE 1 transited 
geosynchronous altitude during a period when high energy 
electrons fluxes were observed. ISEE 1 moved through the 
plasmasheet and was observed to charge to at least -70 volts 
in sunlight. The time of maximum charging, near -100 volts, 
was at 0300 LT . Another geosynchronous satellite, ATS 5, 
was orbiting through the same plasma environment at local 
midnight, and charged to -6 kV in eclipse. The indication 
of charging on ATS 5, indicates that the plasma environment 
was sufficiently hot for spacecraft charging [Ref 13, p. 

5568 to 5578] . Recall that ATS 5 and ATS 6 are two vastly 
different satellites, but charged to similar potentials in 
eclipse (Figure 20) . 

Instrumentation on-board ISEE 1 included two 
electron detectors, the University of Iowa LEPEDEA, and the 
Medium Energy Particle Instrument, or MEPI. The LEPEDEA 
detector data was taken from the radial direction, a 6°x 
38°solid angle. The telemetry for this data covered the 200 
eV to 45 KeV energy range in 32 steps [Frank et al . , 1978]. 
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The MEPI covered the 20 to 1200 KeV energy range 
with eight channels. The detectors provided an angular 
resolution of 10°x 45° with a geometric factor of 10' 2 cm 2 
ster [Williams et al . , 1978]. 

Figure 25 shows the distribution function recorded 
by I SEE 1 on day 76. There are three distinct segments to 
the distribution function that were observed. The low 
energy segment of electron distribution function was 
measured from a few hundred eV to one keV. The average 
temperature was 314 eV and the plasma density was 0.25 cm' 3 . 
Figure 25 shows the electron distribution function. 
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The second electron population can also be described 
as a Maxwellian distribution. The average temperature of 
this group of electrons is 7.7 keV and the density is 0.15 
cm' 3 . The data was collected over an energy range of 4 to 
36 keV. There is a boundary between the two Maxwellian 
distribution functions at 1.7 keV. The values recorded for 
the energy spectrum below 40 keV are typical of the plasma 
found in the plasmasheet. This data is plotted as a series 
of and labelled "low" and "mid" in Figure 26. 
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Above 40 keV, the MEPI data is plotted with "X", and 
is labelled "energetic". The energy spectrum follows an 
inverse power law distribution and is approximately 
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proportional to an E~ 4 distribution function. To smoothly 
connect the mid-Maxwellian and the energetic distribution 
functions, an energetic density of 0.002 cm -3 was fitted. 
This fit agrees with the previous work done [Ref 23]. 

2 . Secondary Electron Differential Energy Flux 
The secondary electron flux produced by this 
distribution function is shown as a function of incident 
energy in Figure 27. 



> 

<v 



V) 

V) 

rst 

£ 

U 

\ 

<L' 



x 

Lx. 

X 

CT> 

<I> 

C 

LxJ 

O 

c 

a> 

<v 



Q 

CT 

O 




Energy (KeV) 



Figure 27 



The majority of the electrons emitted result from the 
differential energy flux below approximately 2 keV. This is 
because the peak of the secondary yield curve is between 300 
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and 400 eV. This corresponds to the peak in the low energy 
segment of the incident flux (Figure 26) . The flux due to 
secondary emission is higher than the ambient flux. For 0 
to 1.5 keV, the peak incident flux is 1 x 10® e/ (cm 2 s sr 
keV) and the peak outgoing secondary flux is 4 x 10® e/ (cm 2 
s sr keV) . 

This reflects the fact that the low energy segment 
of the differential energy flux is between the first and 
second crossover energies for this satellite material, as 
seen in Figure 17. This shows that for low energy electrons 
in the ambient differential energy flux incident onto the 
spacecraft surface, the secondary yield differential energy 
flux over-balances the incident flux. 

3. Backscattered Differential Energy Flux 

The backscattered differential energy flux is 
directly dependent on the incident ambient differential 
energy flux. The particles that are emitted from the 
satellite surface below 50 eV are included in the 
calculation of the secondary yield differential energy flux. 
The two Maxwellian components to the distribution function 
both contribute to backscattered electron flux. Figure 28 
shows the backscattered differential energy flux, as a 
function of energy, for the ISEE 1 spectrum. The 
backscatter flux is an order of magnitude below the ambient 
differential energy flux. The peaks in the backscatter 
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Log Differential Energy Flux (e/crrTs sr KeV) 



curve correspond with the peaks in the incident spectrum 
because the yield curve is relative flat. The net effect of 
backscattering is to reduce the net current due to the 
ambient (0 to 40 keV) electrons by a factor of 10% to 20%. 
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Figure 28 



For more energetic electrons, the influence of 
surface atom interaction will decrease and they will be more 
likely to bury themselves into dielectric surfaces, as shown 
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in Figure 22. The remaining calculations , therefore, are 
associated with correctly determining the deposited charge. 
To do this we need to determine the radiation-induced 
conductivity . 

4 . Radiation-Induced Conductivity 

Incident solar radiation onto dielectric materials 
on spacecraft surfaces will,, over time, cause the RIC of the 
surface to change. The change in conductivity directly 
effects the rate at which charge that enters the dielectric 
material when energetic electrons bury themselves, leaks 
out. This, in turn, helps determine whether enough charge 
can build up in the dielectric to surpass the dielectric 
breakdown threshold and cause an electrical discharge. The 
energetic portion of the differential energy flux is shown 
in Figure 29, in units of e/ (cm 2 s MeV), from 10 to 200 keV. 

Calculation of RIC requires a calculation of dose. 
This is done by convolving the energetic current density, 

J (E) , the collisional stopping power (Figure 23), and the 
specific thickness curve (Figure 24) to obtain the 
differential dose rate (gm sec)" 1 . Integrating the 
differential dose rate from 10 to 450 keV gives the dose 
rate in MeV/ (gm sec) . Figure 30 shows a graph of the 
penetration depth, in cm, of energetic electrons as a 
function of the incident dose rate. 
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Figure 30 
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The figure shows an overlay of different energies 
and how far into the dielectric material they penetrate. 

The lowest energy was 10 keV and the highest was 450 keV. 

The dielectric material used for this graph was borosilicate 
glass. An upper bound of 230 keV was used as it is the 
maximum energy incident onto a .03 cm solar array cover cell 
that will deposit itself into the material. Higher energy 
electrons will pass through the cover cells and into the 
body of the spacecraft. 

The integrated dose rate for the 0.03 cm dielectric, 
is - 0.06 rad/sec. Multiplying this dose rate by the 
largest value in table 5 (k wa2 = 10" 16 sec/ (£2 cm rad) gives 
0 RJtc = 6 x 10" 18 (f2 cm)" 1 . The value of <J L for Corning 7 940 

glass is 1 x 10" 15 (f2 cm)" 1 [Ref 41, p. 141]. The 

conductivity change from radiation is ~ two orders of 
magnitude below the intrinsic conductivity of the 
dielectric. This shows that RIC is an minor factor in the 
overall deposition of charge and, for the purposes of this 
thesis, RIC can be ignored. 

5 . Surface Charge 

The distribution function for day 76, as shown in 
Figure 25, was used to calculate the flux components through 
the satellite surface. Table 6 shows the flux components 
and their respective contributions to the net flux. The 
table is organized as follows. The ambient flux is divided 
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into three components, the low, mid, and energetic energies. 
Summing each contribution gives the incident flux. That 
flux which is produced from secondary emission is labelled 
secondary flux. The backscatter flux is made up of two 
components, the low and mid energy fluxes. The net flux is 
the incident flux minus the sum of the secondary and 
backscattered fluxes. 



TABLE 6 



COMPONENT 


e/m 2 s x 10 11 


Flux < 1.7 keV 


7.24 


Flux > 1.7 keV & < 40 keV 


20.76 


Flux > 40 keV 


1.39 


Total Incident Flux 


29.4 


Secondary Flux 


15 . 9 


BS Flux > .05 keV & <1.7 keV 


. 65 


BS Flux > 1.7 keV & < 40 keV 


2.53 


Total BS Flux 


3.18 


Net Flux 


10 . 3 



The ambient current is .47 (JLA/m 2 . The energetic 
flux alone contributes a current density of 0.022 |!A/m 2 , a 
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very small part of the total. Combining the secondary and 
backscatter fluxes gives the amount of flux that leaves the 
surface. The current density J leaving the surface is 
^backscatter + secondary^ = -306 |iA/m 2 . The net current due to the 
ambient electrons is approximately 0.16 |iA/m 2 negative. 

Hence, a glass surface in shadow should charge negatively. 
The photoelectric saturation current density from Table 4 is 
nominally 21 (IA/m 2 . This shows that when the satellite is 
in sunlight, the current density from incident photons is 
two orders of magnitude higher than the current due to the 
plasma. Therefore, whenever the satellite is in sunlight, 
the amount of charge that leaves the surface is greater than 
that deposited, and the potential of the vehicle should be 
positive . 

If we make the assumption that the satellite exits 
the midnight to dawn sector with a vehicle potential close 
to zero, then the potential will become positive from 
photoelectric emission in the absence of a warm/hot plasma. 
The I SEE 1 period of rotation is three seconds. This means 
that while in sunlight, approximately 20% to 25% of the 
vehicle is in sunlight. Using equation (24) and a spin 
period of three seconds, a net charge of 12.7 [XC/m 2 leaves 
the spacecraft surface each period. This effect dominates 
the other mechanisms and is responsible for maintaining the 
potential close to zero in the absence of an energetic 
plasma . 
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6 . Deep Dielectric Charge 

There is an additional aspect of the charging 
process that must be considered. Energetic electrons of a 
few keV and above began to bury themselves into the body of 
the dielectric, where they become trapped. This buildup of 
charge within the dielectric can develop over an extended 
period of time, and the associated electric fields may 
surpass the breakdown dielectric strength, typically on the 
order of 10 6 V/cm [Ref 34, p. 569 to 591] . For a 0.03 cm 
cover cell, the potential threshold can be as low as 3000 
volts negative, but is more likely to be higher. This has 
been the focus of previous studies of deep dielectric 
charging such as the one by Reagan [Ref 21] . This charge 
deposition has another consequence, however. 

Consider the time for the spacecraft to move through 
the eclipse to be one hour. The deposited charge which 
results from the "energetic" portion of the spectrum depends 
on the lower limit in the integration of J (E) over E. 
Previous work at the Naval Postgraduate School focused on 
Emin = 100 keV. If the lower limit is reduced, a wider 
energy range is included in the integration. This 
substantially increases the energetic current, and hence the 
deposited charge. 

Figure 31 shows J (E) integrated over energy from E min 
to infinity. The data used for the calculation are from 
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Figure 26. Again, ISEE 1 was passing through geosynchr onous 
altitude at 0300 LT [Ref 13] . By varying the lower limit of 
integration, the current density inversely increases. The 
vertical axis is the result of q x I J (E) dE . As you might 
expect, the bulk of the current contribution comes at the 
low end of the energy spectrum. 
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The question, then, is how much charge can be 
deposited over the period of interest. Figure 31 determines 
the time period for the ISEE 1 satellite to charge via our 
proposed mechanism. Conversely, taking one hour as the 
period, then the charge is proportional to the flux and will 
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range from - 180 |iC/m 2 down to - 20 |lC/m 2 . These numbers 
are comparable to, or larger, than the net charge which is 
emitted as a result of photoemission in a spin period. This 
differs substantially from the previous result [Ref 23] 
where the low energy limit was 100 keV. This shows that the 
recommendation for further calculations at lower energies 
was correct. Note that the energy range (10 to 30 keV) not 
included here, will also contribute, though surface 
conductivity will render these less deeply buried electrons 
ineffective for long time scale processes. The inequality 
in equation (24) is satisfied, indicating that the mechanism 
proposed for the development of differential charge on ISEE 
1 satellite should be viable. 

The ISEE 1 charging event was nearly unique for this 
vehicle. The reason appears to be that the energetic 
electron flux was relatively high on this day. Data from a 
second satellite provides selected examples of events with 
very high fluxes of energetic electrons. 

B . SCATHA 

Reagan et al . , surveyed the near— geosynchronous SCATHA 
data set and identified the peak energetic flux events [Ref 
21] . Figure 32 shows all of these high flux days. One of 
these large events was on 29 August, 1979. The data was 
recorded by the P78-2 SCATHA satellite, which was at 6.5 R E , 
in sunlight, and at 0230 LT . The differential flux, as 
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shown in Figure 32 for day 241a, was used here. This 
energetic plasma was then compared to day 76. The energetic 
current density integrated from 40 keV to 1 MeV gives 
^energetic = 0.235 (iA/m 2 . Figure 33 shows the calculations for 
J.n.rg.tic- The energetic plasma of day 241a provided 5 times 
more current than observed on ISEE 1 . Such fluxes would 
cause substantial charge deposition in 5 to 10 minutes. 

This implies that although high, the ISEE 1 fluxes were 
not record setting, and that deep dielectric charging could 
become important on a time scale of ~ 10 minutes. Note that 
on day 241a, the electron flux is very high from 4 to 40 
keV. This energy range is responsible for both deep 
dielectric charging and subsequent surface charging in the 
SCATHA survey. Figure 34 shows the SCATHA potential on day 
241a plotted against time. This event was one of largest 
recorded on SCATHA. 
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Figure 33 
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V 



CONCLUSIONS 



Spacecraft at geosynchronous altitude encounter 
energetic plasmas as they move through the plasmasheet. 

These hot plasmas can effect satellite functions. EMPs 
induced in spacecraft electrical systems can cause data 
upsets, system resets, and even system failures. 

Traditional charging studies focused on the nearly 
instantaneous processes found in eclipse [Ref 7], and 
daylight charging processes with time constants of minutes 
[Ref 17, Ref 8]. This work has shown that deep dielectric 
charging, which occurs on time scales of many minutes, may 
also be significant for daylight charging. 

Deep dielectric charging, occurs when the distribution 
function of the hot plasma contains electrons whose average 
temperature exceeds about 10 keV. This mechanism deposits 
energetic electrons deep within the dielectric materials on 
the satellite surface. The charge deposited within the 
dielectric remains stationary, leaking out over a period of 
many hours. This buildup of charge within the dielectrics 
is independent of the other processes involving the surface. 

The critical question is, if more charge can be 
deposited in an hour or less than emitted from the surface 
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in one spin of the satellite, will an isolated segment of 
dielectric develop a net negative charge? 

Records from very energetic plasmas were examined. Data 
from ISEE 1 daylight negative potential event and the record 
daylight event for SCATHA were studied. The ISEE data 
showed that if 40 to 250 keV electrons are considered, 
sufficient charge is deposited in an hour to result in a net 
negative charge. The five times higher fluxes observed on 
SCATHA would reduce the charging time scale to a few 
minutes . 

Deep dielectric charging is important and may determine 
the daylight potential behavior of geosynchronous 
satellites . In particular, it may provide a charging 
mechanism for rapidly spinning satellites, which do not 
otherwise have the tendency to differentially charge. 
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VI . RECOMMENDATIONS 



Follow up work should focus on the mechanism by which 10 
keV to 200 keV electrons bury themselves into dielectric 
materials. Modelling of the process would produce numerical 
values that should be compared to recorded charging events 
and the attendant electron distribution function. The 
speculation that deep dielectric charging will lead to 
daylight negative potentials needs to be studied via 
modelling . 
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